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Abstract-  The  Poincare  plot  is  an  emerging  Heart  Rate  Vari¬ 
ability  (HRV)  analysis  technique,  the  geometry  of  which  has 
been  shown  to  distinguish  between  healthy  and  unhealthy  sub¬ 
jects  in  clinical  settings.  The  Poincare  plot  is  able  to  display 
nonlinear  aspects  of  the  interval  sequence  and  is  therefore  of 
interest  in  characterizing  the  nonlinear  aspects  of  HRV.  The 
problem  is,  how  do  we  quantitatively  characterize  the  geometry 
of  the  plot  to  capture  useful  descriptors  that  are  independent  of 
existing  HRV  measures?  In  this  paper,  we  investigate  a  popular 
existing  category  of  techniques  and  show  that  they  measure  lin¬ 
ear  aspects  of  the  intervals  which  existing  HRV  indices  already 
specify.  The  fact  that  these  methods  appear  insensitive  to  the 
nonlinear  characteristics  of  the  intervals  is  an  important  finding 
because  the  Poincare  plot  is  primarily  a  nonlinear  technique. 

Keywords-.  HRV,  Poincare  plot,  nonlinear  analysis 

I.  Introduction 

The  field  of  heart  rate  variability  (HRV)  studies  the  fluc¬ 
tuations  in  the  intervals  between  heartbeats,  known  as  RR 
intervals.  The  Poincare  plot,  a  technique  taken  from  nonlin¬ 
ear  dynamics,  portrays  the  nature  of  these  fluctuations 
graphically.  It  is  a  scatter-plot  of  each  RR  interval  plotted 
against  the  next  interval.  Poincare  plot  analysis  is  an  emerg¬ 
ing  quantitative-visual  technique  whereby  the  shape  of  the 
plot  is  categorized  into  functional  classes  that  indicate  the 
degree  of  heart  failure  in  a  subject  [1,  2].  The  plot  provides 
summary  information  as  well  as  detailed  beat-to-beat  infor¬ 
mation  on  the  behavior  of  the  heart  [3]. 

Support  is  increasing  for  nonlinear  analysis  techniques  and 
quantitative  descriptors  as  it  has  become  evident  that  the  car¬ 
diac  systems  are  nonlinear  in  their  function  [4].  The  Poincare 
plot  is  becoming  a  popular  technique  due  to  its  simple  visual 
interpretation  and  its  proven  clinical  ability  as  a  predictor  of 
disease  and  cardiac  dysfunction  [5].  The  problem  regarding 
Poincare  plot  use  has  been  the  lack  of  obvious  quantitative 
measures  that  characterize  the  salient  features  of  the  plot. 
Researchers  have  put  forward  a  number  of  techniques  that 
attempt  to  quantitatively  summarize  the  plot’s  geometric  ap¬ 
pearance.  The  efforts  can  be  summarized  into  3  categories: 
geometrical  descriptors,  scanning  parameters  and  image  dis¬ 
tribution  parameters  [6].  Of  these,  the  geometrical  descrip¬ 
tors  are  the  most  popular  in  the  clinical  and  physiological 
HRV  literature. 

In  this  study,  we  consider  the  geometrical  Poincare  plot 
descriptors.  We  provide  expressions  that  connect  each  de¬ 
scriptor  to  existing  linear  measures  of  HRV.  This  accom¬ 
plishes  two  things.  Firstly,  it  provides  insight  into  Poincare 
plot  geometry  in  terms  of  the  well-understood  existing  indi¬ 
ces  of  HRV.  Secondly,  it  shows  that  these  measures  are  not 
independent  to  the  existing  standard  linear  statistics.  There¬ 
fore,  the  intrinsic  ability  of  the  Poincare  plot  to  identify  non¬ 


linear  beat-to-beat  structure  is  not  being  exploited  by  these 
techniques. 

II.  The  linear  HRV  indices 

This  section  describes  the  standard  linear  indices  of  HRV. 
In  this  paper,  the  time-course  of  the  RR  intervals  is  denoted 
by  RR„,  with  n  —  \..N  .  We  assume  wide-sense  stationarity 
so  the  following  basic  properties  hold:  E[RRn  ]  =  E[RRn+m  ] 
and  E[RR;]  =  E[RRl J. 

A.  Standard  deviation  of  the  RR  intervals 

The  standard  deviation  of  the  RR  intervals,  denoted  by 
SDRR  ,  is  often  employed  as  a  measure  of  overall  HRV.  It  is 
defined  as  the  square  root  of  the  variance  of  the  RR  intervals, 

SDRR  =  (1) 

where  the  mean  RR  interval  is  denoted  by  RR  =  E  [RR„  ] . 

B.  Standard  deviation  of  the  successive  differences 

The  standard  deviation  of  the  successive  differences  of  the 
RR  intervals,  denoted  by  SDSD  ,  is  an  important  measure  of 
short-term  HRV.  It  is  defined  as  the  square  root  of  the  vari¬ 
ance  of  the  sequence  A RRn  =  RRn  -  RRn+l  (the  delta-RR  in¬ 
tervals), 

SDSD  =  ^E[ARR;]  - ARRn2  .  (2) 

Note  that  ARRl:  =  E[RRn]~ E[RRtl+]  ]  =  0  for  stationary 
intervals.  This  means  that  SDSD  is  equivalent  to  the  root- 
mean-square  of  the  successive  differences,  denoted  RMSSD  . 

C.  Autocorrelation  and  autocovariance 

The  autocorrelation  function  is  an  important  measure  of 
HRV  simply  because  its  Fourier  transform  is  the  power 
spectrum  of  intervals.  The  autocorrelation  function  of  the  RR 
intervals  is  defined  as, 

yRR(m)  =  E[RRnRRnlm).  (3) 

Spectral  analysis  is  normally  performed  on  the  mean- 
removed  RR  intervals,  and  therefore  the  mean-removed  auto¬ 
correlation  function,  called  the  autocovariance  function,  is 
often  preferred: 

0™  On)  =  E[(RRn  -  RR)(RR,1+m  -RR)].  (4) 
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The  autocovariance  function  is  related  to  the  autocorrela¬ 
tion  function  for  stationary  intervals  by  the  relationship 

(pRR  (m)  —  yRR  (m)  —  RR  .  The  variance  of  the  RR  intervals  is 
SDRR 2  =<f)RR( 0) ,  and  the  variance  of  the  delta-RR  intervals 
is  SDSD2  =2(fi>RR(fi)-(f)RR(l)) .  Accordingly,  these  indices 
are  linear  measures  of  HRV. 

111.  The  Geometrical  Poincare  plot  descriptors 

The  RR  interval  Poincare  plot  typically  appears  as  an 
elongated  cloud  of  points  oriented  along  the  line-of-identity 
(see  Fig.  1).  The  dispersion  of  points  perpendicular  to  the 
line-of-identity  reflects  the  level  of  short-term  variability  [5]. 
The  dispersion  of  points  along  the  line-of-identity  is  thought 
to  indicate  the  level  of  long-term  variability.  The  geometric 
descriptors  use  statistical  moments,  usually  directed  along  the 
line  of  identity,  to  measure  the  shape  of  the  plot.  In  this  sec¬ 
tion,  we  relate  several  geometrical  techniques  that  are  popu¬ 
lar  in  the  literature  to  linear  measures  of  HRV. 


RR(n+1) 


A.  Ellipse  fitting  technique 

By  far  the  most  popular  technique  used  to  characterize  the 
shape  of  the  plot  numerically  is  the  technique  of  fitting  an 
ellipse  to  the  plot  [6-9],  as  Fig.  1  details.  A  set  of  axes  ori¬ 
ented  with  the  line-of-identity  is  defined  [9],  The  axes  of  the 
Poincare  plot  are  related  to  the  new  set  of  axes  by  a  rotation 
of  9  =  n ! 4  radians: 


-sin# 

cos  0 

In  the  reference  system  of  the  new  axes,  the  dispersion  of 
the  points  around  the  x,  axis  is  measured  by  the  standard 
deviation  denoted  by  SDl  [9].  This  quantity  measures  the 
width  of  the  Poincare  cloud,  and  therefore  indicates  the  level 
of  short-term  HRV  [2,  3,  5,  9].  The  length  of  the  cloud  along 
the  line-of-identity  measures  the  long-term  HRV  and  is 
measured  by  SD2  which  is  the  standard  deviation  around  the 
x2  axis  [2,  3,  5,  9].  These  measures  are  related  to  the  stan¬ 
dard  HRV  measures  in  the  following  manner: 

SOI2  =  Var (x, )  =  Var (j, RRn  -^RR„+l) 

,  (6) 

=  j  Var  (RR„  -  RRII+I )  =  j  SDSD2 

Thus,  the  SDl  measure  of  Poincare  width  is  equivalent  to 
the  standard  deviation  of  the  successive  intervals,  except  that 
it  is  scaled  by  1  /  s[2  .  This  means  that  we  can  relate  SDl  to 
the  autocovariance  function: 

SD12=0jm(  O)-0m(1).  (7) 

It  may  also  be  shown  that  the  length  of  the  Poincare  cloud 
is  related  to  the  autocovariance  function: 


Fig.  1:  An  example  Poincare  plot,  also  detailing  the  ellipse 

fitting  process  and  the  histograms  derived  from  the  plot. 

SD22  =  (fig#  (0)  +  (jtgg  (1)  .  (8) 

By  adding  equations  (7)  and  (8)  together,  we  obtain  the  re¬ 
sult 

SDl2  +SD22  =  2(/>RR  (0)  =  2 SDRR2 .  (9) 

Finally, 

SD2 2  =  2  SDRR2  -  j  SDSD 2 .  (10) 

Equation  (10)  allows  us  to  interpret  SD2  in  terms  of  ex¬ 
isting  indices  of  HRV  and  can  also  be  used  to  argue  that 
SD2  reflects  the  long-term  HRV,  however,  we  delay  the 
discussion  until  later.  Clearly,  fitting  an  ellipse  to  the  Poin¬ 
care  plot  does  not  generate  indices  that  are  independent  of  the 
standard  time  domain  HRV  indices.  In  fact,  the  width  of  the 
Poincare  plot  is  a  linear  scaling  of  the  most  common  statistic 
used  to  measure  short-term  HRV,  the  SDSD  index.  In  other 
words,  the  width  of  the  Poincare  plot  should  correlate  ex¬ 
tremely  highly  with  other  measures  of  short-term  HRV,  as 
indeed  it  does  [3]. 

B.  Histogram  techniques 

Another  method  to  quantify  the  shape  of  the  Poincare  plot 
is  to  measure  the  statistical  properties  of  various  projections 
of  the  plot  via  histogram  distributions  [2,  3,  5].  Fig.  1  shows 
the  three  main  projections  used.  They  are: 

RR  interval  histogram.  The  histogram  of  the  Poincare  plot 
points  projected  onto  the  x-axis  (or  the  y-axis).  This  histo¬ 
gram  is  usually  quantified  by  the  mean  and  standard  devia¬ 
tion,  which  correspond  directly  to  the  standard  linear  meas¬ 
ures  RR  and  SDRR  .  This  view  provides  summary  infor¬ 
mation  on  the  overall  HRV  characteristics. 
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‘Width’,  or  delta-RR  interval  histogram.  The  Poincare 
plot  points  are  projected  along  the  direction  of  the  line-of- 
identity.  ft  is  not  exactly  equivalent  to  the  delta-RR  interval 
histogram  as  the  abscissa  has  been  scaled  by  the  factor 
1  /  V2  .  Mathematically,  it  is  the  distribution  of  x, .  There¬ 
fore,  the  standard  deviation  of  the  width  histogram  is  equal  to 
SDl .  This  histogram  provides  summary  information  on  the 
short-term  characteristics. 

‘Length  ’  histogram.  The  Poincare  plot  points  are  projected 
onto  the  line-of-identity.  The  histogram  is  described  mathe¬ 
matically  by  the  distribution  of  x2  and  the  standard  deviation 
is  therefore  equivalent  to  SD2  .  Consequently,  the  length 
histogram  portrays  the  long-term  characteristics  of  HRV. 
The  dispersion  properties  of  these  histograms  are  character¬ 
ized  by  SDRR  ,  SDl  and  SDl .  Hence,  they  are  linked  to 
the  standard  linear  time-domain  measures  of  HRV 

C.  Correlation  coefficient 

Some  researchers  have  employed  the  correlation  coeffi¬ 
cient  of  the  Poincare  plot  to  characterize  its  shape  [10].  This 
measure  is: 


E 

\RR„-RR)(RRn+l-RR) 

£ 

(RR„-RR)2^E 

\(RRn+1-RR)2] 

(ID 


For  the  Poincare  plot,  the  correlation  coefficient  can  be  ex¬ 
pressed  in  terms  of  the  autocovariance  function: 


RR(n+l) 

b 


•  (b.a) 


a)  b) 

Fig.  2:  a)  Poincare  plot  of  alternating  sequence,  b)  Diagram 
of  a  third-order  Poincare  plot. 


SDRR  is  not  zero.  We  propose  that  Poincare  plot  length  is  a 
more  consistent  and  appealing  measure  of  long-term  vari¬ 
ability  as  the  length  of  the  Poincare  plot  is  zero  for  an  alter¬ 
nating  sequence  indicating  zero  long-term  variability,  as  Fig. 
2a  shows. 

With  these  ideas  in  hindsight,  we  can  now  explain  the  sig¬ 
nificance  of  equations  (9)  and  (10).  Equation  (9)  states  that 
the  sum  of  the  short-term  and  the  long-term  variability  is  the 
total  variability.  Equation  (10)  affirms  that  the  long-term 
variability  is  the  total  variability  minus  the  contribution  due 
to  short-term  variability. 


IV.  Generalizations  of  the  Poincare  plot 

Two  different  types  of  scatter-plots  encountered  in  the  lit¬ 
erature  can  be  considered  simple  generalizations  of  the  Poin¬ 
care  plot. 


o»=<MD/<M0).  (12) 

Therefore,  the  correlation  coefficient  is  a  linear  measure, 
even  though  it  is  based  on  the  Poincare  plot  which  displays 
nonlinear  features.  None  of  the  geometric  descriptors  that 
have  been  analyzed  in  this  section  are  sensitive  to  the  nonlin¬ 
ear  features  the  Poincare  plot  displays. 

D.  Short-  and  long-term  variability 

The  relationships  identified  so  far  give  clear  mathematical 
insight  into  what  the  length  and  width  of  the  Poincare  plot 
depict.  The  length  and  width  of  the  Poincare  plot  have  been 
suggested  as  indicative  of  the  levels  of  long-  and  short-term 
variability.  It  is  reasonably  clear  that  the  standard  deviation 
of  the  delta-RR  intervals,  as  measured  by  SDSD  ,  RMSSD 
or  SDl ,  is  a  measure  of  short-term  HRV.  In  fact,  this  state¬ 
ment  can  be  made  even  more  precise:  these  indices  are  meas¬ 
ures  of  the  variability  over  a  single  beat.  The  standard  devia¬ 
tion  of  the  RR  intervals,  as  measured  by  SDRR  ,  is  often  em¬ 
ployed  as  a  measure  of  long-term  HRV.  However,  there  is  a 
problem  with  this  interpretation  because  this  quantity  meas¬ 
ures  all  the  variability,  long-term  and  short-term. 

For  example,  take  a  set  of  RR  intervals  that  has  variability 
only  over  a  single  beat.  Such  a  sequence  can  be  described  as 
alternating  between  two  values,  e.g.,  a,  b,  a,  b,  etc.  It  is  clear 
that  the  sequence  contains  variability  only  over  a  single  beat: 
every  second  beat  is  equivalent.  If  SDRR  is  taken  to  be  an 
index  of  long-term  variability,  a  contradiction  arises  because 


A.  Lagged  Poincare  plots 

Instead  of  plotting  RRn  against  RRn+l ,  some  researchers 
have  investigated  plotting  RRn  against  RR„+m  where  m  is 
allowed  to  vary  from  1  to  some  small  positive  value.  In  gen¬ 
eral,  the  plot  is  still  clustered  around  the  line-of-identity. 
However,  the  length  and  width  of  the  plot  are  altered  as  the 
lag  is  increased.  It  is  possible  to  show  the  width  and  length 
measures  SDl  and  SDl  can  be  generalized  for  lag  m  : 

SDl(mf  =<pRR(0)-<pRR{m) 

SDlfmf  =<l>RR(0)+<l>RR(m) 

The  length  and  width  of  the  lag-  m  Poincare  plot  is  related 
to  the  covariance  function  at  lag  m  .  Note  also  that: 

<pRR(m)  =  \(SDl(m)2 -SDlfnf)  (14) 

This  result  is  very  interesting,  as  it  states  that  the  set  of 
lagged  Poincare  plot  length  and  widths  are  a  complete  de¬ 
scription  of  the  autocovariance  function,  and  hence,  also  the 
power  spectrum  of  the  intervals.  Equation  (14)  also  provides 
us  with  a  geometrical  relationship  between  the  autocovari¬ 
ance  function  and  the  Poincare  plot’s  shape.  If  <fRR(m)  =  0 
then  SDl  =  SDl  and  the  length  and  width  of  the  plot  are 
equal.  If  (j)RR(m)  >  0  then  SDl  <  SDl  and  the  plot  is  longer 
than  it  is  wide,  i.e.  dominated  by  short-term  activity  and  vice 
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versa  for  long-term  variability.  This  is  similar  to  the  concept 
of  a  scatter  plot’s  correlation  coefficient  sign.  In  fact  the 
series  of  correlation  coefficients  of  the  lagged  Poincare  plots 
are  simply  a  scaled  version  of  the  autocovariance  function: 

rRRi<m)  =  <pRR(m)l<pRR(  0)  (15) 

B.  Higher-order  Poincare  plots 

The  standard  Poincare  plot  is  a  scatter-plot  of  the  pairs 
(RRn ,  RRn+1 ) ,  and  is  considered  to  be  of  first  order.  The  sec¬ 
ond  order  Poincare  plot  is  a  three-dimensional  scatter-plot  of 
the  triples  (RRn,RRn+l,RRn+2) .  There  are  three  orthogonal 
views  of  the  three-dimensional  shape  of  this  plot,  each  being 
a  view  along  one  of  the  axis.  These  views  result  in  two- 
dimensional  projections  of  the  three-dimensional  cloud  onto 
each  of  the  coordinate  planes  ( RRn , RRn , , ) ,  (RRlt+l,RRlt+2) 
and  (RRn,RRll+2) .  The  first  two  views  are  equivalent  to  the 
standard  Poincare  plot  and  the  third  is  the  lag-2  Poincare  plot 
as  Fig.  2b  shows.  This  idea  can  be  extended  into  higher  di¬ 
mensions,  with  the  projections  of  the  plot  onto  coordinate 
planes  being  lagged  Poincare  plots.  So,  an  order  m  Poincare 
plot  is  geometrically  described  by  the  set  of  lagged  Poincare 
plots  up  to  and  including  lag  m  .  The  results  for  lagged 
Poincare  plots  carry  over  to  higher-order  Poincare  plots. 

V.  Conclusion 

We  have  shown  that  the  so-called  geometric  techniques 
that  characterize  the  geometry  of  a  Poincare  plot  are  related 
to  linear  indices  of  HRV.  In  addition,  we  provide  arguments 
supporting  the  claim  that  the  width  of  the  Poincare  plot  corre¬ 
sponds  to  the  level  of  short-term  HRV,  while  the  length  of 
the  plot  corresponds  to  the  level  of  long-term  variability. 

The  methods  of  quantifying  the  Poincare  plot  that  we  have 
investigated  herein  are  not  capable  of  depicting  the  additional 
beat-to-beat  variability  information  shown  on  a  Poincare  plot. 
However,  the  additional  information  is  likely  to  be  of  consid¬ 
erable  value.  The  fact  that  that  the  length  and  width  of  a 
Poincare  plot  corresponds  so  conveniently  to  standard  time 
domain  statistics  is  a  very  nice  feature.  However,  simply 
treating  the  Poincare  plot  as  a  tool  for  graphically  represent¬ 
ing  these  statistics  is  to  ignore  some  of  its  potential  capabili¬ 
ties.  Therefore,  we  recommend  nongeometric  techniques, 
such  as  scanning  parameters  [7,  11]  and  image  distribution 
measures  [12],  be  investigated  for  a  full  appraisal  of  the 
clinical  capabilities  of  the  Poincare  plot. 
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